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Abstract 

We consider a single disk moving under the influence of a 2D viscous fluid and we study the asymptotic 
as the size of the solid tends to zero. 

If the density of the solid is independent of e, the energy equality is not sufficient to obtain a uniform 
estimate for the solid velocity. This will be achieved thanks to the optimal L v — L q decay estimates of the 
semigroup associated to the fluid-rigid body system and to a fixed point argument. Next, we will deduce 
the convergence to the solution of the Navier-Stokes equations in R 2 . 


1 Introduction 


We study in this paper the asymptotic of a fluid-solid system as the solid is a rigid disk which shrinks to a 
point. We first describe the fluid-solid system. 

1.1 The fluid-solid system 

We consider a rigid disk 


S e (t) = B(h s (t),e) 

immersed into a viscous incompressible fluid. At time t > 0, the lagrangian coordinates to the body read 

v(t, x) := h e (t) + R gs(t) (x - h 0 ), 

where for all t, 

h E {t) e R 2 with h e ( o) = h 0 , e e (t) ei, R e = ( cos6 q ~ sm ^ . 

w w w ysm 0 cost/ J 

The domain of the fluid evolves through the formula 

R E (t) :=R 2 \S E (t). (1.1) 

We denote by n := n e (t,x ) the exterior unit normal of d!F E (t). The equations for the fluid-solid system read 


<9?/ £ 

—-1- (u e ■ X7)u e — div cr(u e ,p e ) = 0 

at 

t > 0, x e R E (t), 

(1.2) 

div u e = 0 

t > 0, x e R E (t), 

(1.3) 

O 

II 

-f-T 

to 

S3 

t > 0, 

(1.4) 

| X | — ^oo 



u e = (h e )'(t) + {6 E )'(t)(x-h e (t)) A - 

t > 0, x £ dS E (t), 

(1.5) 

m e {h e )"{t) = — f a(u e ,p e )n dj 

t > 0, 

(1.6) 

J dS e (t) 



J e (9 e )"(t) = — f (x — h e ) ± ■ a(u E ,p e )n dj 

t > 0, 

(1.7) 

J dS s (t) 



u e { 0,-) = ul 

in T E , 

(1.8) 

h E (0) = h o , (h E )'(Q)=ll, 9 E { 0) = 0, (0 e ) , (O)=r§. 

(1.9) 
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Here and in what follows 

with v > 0 is the constant viscosity and 


cr(u,p) = 2 vD(u) — pi 2 , 


Din) := \ ((Vu) + (Vu)*). 


We write for any x G 


x 1 - := 


-x 2 

Xl 


It 7r/2*^' 

It is convenient to extend the velocity field u £ inside the rigid disk as follows: 

u £ (t, x) = ( h £ )'(t ) + ( 9 £ )'(t)(x — h £ (t)) ± t > 0, x G S £ (t). 


( 1 . 10 ) 


To apply the result in m , we also need to assume that the center of the mass corresponds with the center 
of the disk. For simplicity, let us assume that the density p e > 0 is constant in the disk. We define a global 
density in R 2 by 


p £ (t,x) 


1 x £ .F £ (f), 
p £ x G S £ (t). 


t > 0. 


For any smooth open set O, we define 


• V(0) := \p G Cg°(0) | divv? = 0 in C>}; 


• H(0) the closure of V(O) in the norm L 2 : 


11(0) = G L 2 (0) | divy: = 0 in O, p ■ n = 0 at <9oj; 


• V(O) the closure of V(0) in the norm H 1 : 

V(0) = {p£ Hq(O) | divv? = 0 in C>} 

and its dual space by V(O) with respect to 11(0). 

We also define V R (F £ (t)) the subspace of V(R 2 ) of velocity fields that are rigid in the solid: 

V R (T £ (t)) := {p G H\R 2 ) ■ D(p) = 0 in S £ (t), div^ = 0} . (1.11) 

Under the following hypotheses on the initial conditions 

Uq G L 2 (J ? o), divrtg = 0, Uq ■ n = 1% ■ n on 6Sq, (1.12) 

there exists a unique global weak solution (u £ , h £ , 9 £ ) see [M] , in the sense of the definition below. 

Definition 1.1. We say that ( u £ , h £ , 0 £ ) is a global weak solution of (TOl) OP if, for any T > 0, we have 

u £ G L°°(0, T; L 2 (R 2 )) n L 2 (0,T; H\M. 2 )), u £ (t , •) G V R (T £ (t)), 
and if it satisfies the weak formulation 

/ f p E u £ ■ (—p— + (u £ • V)</? e ^) dxds + 2v f f D(u £ ) : D(p e ) dxds = f p £ u £ 0 (x) ■ p E (Q , x) dx, (1.13) 

J 0 J R 2 \ ot ) J o J K 2 Jr 2 

for any p £ G C*([0, T); 17 1 (R 2 )) such that p £ (t,-) G V R (iF e (f)). 
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1.2 Massless pointwise particle in the whole plane 

When e —> 0, we establish the convergence of u e to the unique solution of the Navier-Stokes equations in the 
whole plane R 2 . 

The asymptotic behavior of the fluid motion around shrinking obstacles is already considered in several 
recent papers. Iftimie, Lopes Filho and Nussenzveig Lopes (21 have studied the case of one small fixed obstacle 
in an incompressible viscous fluid in 2D. Iftimie and Kelliher [2D] have treated the same situation in 3D. In 
[2T[;25] Lacave has considered the case of one thin obstacle shrinking to a curve in 2D and 3D. 

There is also a large literature about porous medium in the homogenization framework. Since the pioneer 
work of Cioranescu and Murat for the Laplace problem, the Navier-Stokes system was studied, in particular, 
by Allaire mm- We also mention [3 01271 ED! HUES] for the fluid motion through a perforated domain. 

In all the above studies, the general strategy relies on energy estimate to get a uniform estimate in H 1 . 
It turns out that such an estimate is sufficient to pass to the limit in the weak formulation by a troncature 
procedure. Namely, for a test function <p £ T>(Ll) and for a cutoff function % £ , we note that x e T is an admissible 
test function for the Laplace problem in the perforated domain fi e . If the inclusions are far enough, for standard 
cutoff function, ||x e II w 1 -** remains bounded only for p 2 in dimension two, which allows to pass to the limit 
in terms like f Vu e : V(x e p). For the Navier-Stokes equations, the cutoff procedure is more complicated and 
relies on Bogovskii operators. Indeed, we need approximated test functions that are divergence free (see 0 [251 
and Section IQ) . 

When the obstacles can move under the influence of the fluid, we also need to control uniformly the velocities 
of the solids. In the case of the system (03 00) or more generally in the case of a system with several rigid 
bodies moving into a viscous incompressible fluid, this control can be obtained from the energy estimate if the 
masses are independent of e (see Section 15721 for details). 

If the masses tend to zero, it is no more possible to deduce estimates of the velocities of the solids indepen¬ 
dently of e from the energy estimate. One could try to get an estimate of rigid velocities from the boundary 
condition. However, since the size of the solids tend to zero, this leads to look for a C'°-estimate for the fluid 
velocity, and thus for H s estimates with s > 1. It was the strategy followed in mm with a H 2 analysis. 
Unfortunately, these articles are based on uniform elliptic estimates in the exterior of a small obstacle which 
fail for s > 1 (see a counter-example related to these estimates in 0). 

Our strategy is different here. Our basic remark is that the small obstacle problem is related to the long-time 
behavior though the scaling property of the Navier-Stokes equations u E (t,x) = e~ 1 u 1 {e~ 2 t, e~ l x). For one disk 
moving in the plane, the long-time behavior has been recently studied by Ervedoza, Hillairet and Lacave in 
|12j . In particular, they have obtained the optimal decay estimates of the Stokes semigroup, i.e. with the rates 
corresponding to the heat kernel (which are invariant to the parabolic scaling). These estimates are the key to 
treat the massless pointwise particle. 

The goal of the main theorem is to treat the case where the disk shrinks to a massless pointwise particle : 

P e = Po (1-14) 

hence, 

m £ = e 2 m l and J e = e 4 J 1 . (1-15) 

We consider the massless case for small data: 

Theorem 1.2. A ssume (11.141) . Then there exists Ao such that the following holds. 

Let (uo>^o> r o) a family in L 2 (J- jjj) xl 2 xl verifying (11.121) and such that 

e Kol> e 2 kol> II u oIIl 2 (^) < A 0 (1.16) 

and 

Uq —*• uq in L 2 (R 2 ). (1.17) 

Then for any T > 0 we have 

u e ^u in L oo (0,T;L 2 (R 2 ))nL 2 (0,T;i7 1 (R 2 )) (1.18) 
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where u is the weak solution of the Navier-Stokes equations in R 2 associated to uq: for any ip £ Cl ([0, T); V(R 2 )), 

[ u ■ ( + ( u ■ V)<p ] dxds + v f f Vu : "Vipdxds = f uo(x) ■ ip( 0, x) dx. 

R 2 \ at ) J g J R 2 Jr2 

Remark 1.3. For any T > 0, we will actually establish in Section B~T1 that. up to a subsequence, we have 

h E —>• h uniformly in [ 0 ,T], 

and in Section H~3l that, for any Ogl 2 \ {h(t)} (for all t £ (ti,£ 2 ) with some t\,t 2 £ [0,T]), we have 

Pqu e —> P o u strongly in L 2 (t\,t 2 ', L 4 (0)), 

where P q is the Leray projector. This strong limit will be used to pass to the limit in the non-linear term. 

As we recover at the limit the weak solution of the Navier-Stokes equations in the whole plane, and as this 
solution is unique by the Leray theorem, we will deduce that we do not need to extract a subsequence in (11.181) . 

For a 2D ideal incompressible fluid governed by the Euler equations, the case of a massive pointwise particle 
(i.e. where m E = m l is independent of e) in the whole plane was treated in E, a massless pointwise particle in 
the whole plane in nu and both case in a bounded domain in jlSl . In these works, non-trivial limit was obtained 
(namely, Kutta-Joukowski lift force or vortex-wave system) when we consider non-zero initial circulations around 
the small solids. 

1.3 Plan of the paper 

The remainder of this work is organized in four sections. 

In the next section, we provide three examples where the initial convergence (11.171) holds. 

We establish in Section 0 some uniform estimates on ( u e , ( h e )'). The energy estimate will give us directly a 
good estimate for the fluid velocity u e but not for the disk velocity (h E )'. Thanks to the results of [T?, we will 
prove that some L p — L q estimates of the Stokes semigroup are independent of e, and by a fixed point argument 
we will get a uniform estimate of the disk velocity. 

Section [3] is dedicated to the passing to the limit. We introduce the cutoff procedure which follows the 
trajectory of the solid. A crucial point is to construct a corrected test function (p v which satisfies the divergence 
free condition. This will be obtained by the Bogovskii operator El HI- Then we follow the analysis developed in 
[25] . Roughly, we pass first to the limit e —> 0 far away from the solid to get that u satisfies the Navier-Stokes 
equations in this region. Next, we pass to the limit rj 0 in the cutoff function, to prove that the equations 
are also verified in the vicinity of the massless pointwise particle. 

We will discuss in the last section the three dimensional case and we will give the extension of our main result 
in the case where several solids (with any shapes) tend to massive pointwise particles in bounded domains. In 
this case the energy estimate is sufficient to obtain a uniform estimate of the solid velocities, which allows us to 
reach more general geometric configurations than in the massless case. 

2 Examples of initial conditions 

In this paragraph, we develop three examples of family (uq,£q) e satisfying the compatibility condition (11.121) 
which converges in L 2 (R 2 ). 

Example 2.1. The first trivial example is the case where Uq is independent of e. Namely, let us consider 
(uq° , , rg°) satisfying (11.121) in for some £0 > 0. The extension of Uq° by £q° + rQ°(x — ho) 4 - inside the 

disk verifies also ( 11 . 121 ) in Eq for any e £ ( 0 ,£o]. 

Example 2.2. In domains depending on £, a standard setting is to give an initial data in terms of an independent 
vorticity uiq = curlrtQ and initial circulation 70 = §oB(h 0 e) u o ' T ( see i e -S-; PH Q 21 [201 ED [Ml EH)- For the 
2D Euler equations, the vorticity is the natural quantity because it satisfies a transport equation, which implies 
some conserved properties (e.g. the L p norm of the vorticity for p £ [1, 00] and the circulation of the velocity). 
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For the 2D Navier-Stokes equations in domains with fixed boundaries, the vorticity and the circulation are 
less relevant, because the Dirichlet boundary condition implies that the circulation is zero for t > 0, and the 
vorticity equation does not give anymore the conservation of the L p norm. For these equations, the standard 
framework is related to the energy estimate (13.121) . i.e. to consider initial data belonging to L 2 (Tq). In terms 
of the vorticity, we recall that Ug(x) = (70 + hz Ct , o) 2 ^| x p + C>(|^p) at infinity (see for instance [HI Section 2] 
or m Section 3]), hence 

u e 0 G L 2 (Tq) ■<=> 70 + [ wo = 0. 

J *5 

For these reasons, the most natural condition is 70 = / w 0 = 0. In this case, it is easy to prove the following 
result. 

Lemma 2.3. Let (£o,ro) G R 3 , w 0 G L“(R 2 \ {0}) fixed such that f u> 0 = 0. Then, for e small enough such 
that supp wq D B( 0,e) = 0, we have a unique solution Ug in L 2 [T,g) of 


div«Q = 0 in Tg, curing = wo in £Fg, lim Ug(x) = 0, 

| | —>-oo 


Zq • n = £q ■ n on dB{ 0, e), (7) 

fa 


8B( 0,e) 


Uq ■ T ds = 0. 


Moreover, extending Ug by £q + rox 2 - in B( 0,e) we have 


Uq —*• ug weakly in L 2 (R 2 ), 


where ug = A' R 2 [wq] = 2 ^\ x \ 2 * * s vector field in L 2 (M 2 ) such that 


divito = 0 in R 2 , curlrto = wo in R 2 , lim uo(x) = 0. 

| X | —^OO 

Proof. The existence and uniqueness of Ug is well-known (see e.g. m Section 2]): 




B(0,s) c 


(x-y) 1 - , w , 1 

-w 0 (y)dy + 


x x — e 2 y* \ 


'-y I 5 


2-7T 


B( 0,e)° 


\x\ 2 \x — e 2 y*\ 2 


) w 0 (y)dy 


with the notation y* = y/\y\ 2 . By a standard computation, we note that 



£ 


2 



L 2 (^S) 


< Ce \£ 0 |. 


It is also rather classical to prove that the second integral in the right hand side tends to zero as e —> 0. For 
instance, the authors establish in ED Lemmas 7 and 10] the uniform estimate in L 2 (R 2 ) and the convergence 
to zero in V'(R 2 ), which implies the weak convergence in L 2 (R 2 ). 

As h Ib( o,e) c u °(y) d y = = u o, this ends the proof. □ 

Remark 2.4. Even if the zero circulation condition is mandatory for strong solutions to the Navier-Stokes 
equations in fixed domain, for the fluid-solid system the no-slip boundary condition would imply 


i 

J dB{ho,e) 


Uq - t ds = 27 t e 2 rg. 


Hence, an interesting extension could be to study the case of non zero initial circulation 7 ( 5 . If we assume that 
7 q is independent of e, some singular terms appear at the limit of the form 7 o 2 ^-°h 0 \ 2 > which does not belong to 
L 2 oc (R 2 ), but only to Lf oc (R 2 ) for all p G [1,2). Another difficulty in this case is to ensure that 7 g + w 0 = 0 
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for any e, in order to state that the initial velocity is square integrable at infinity. A possibility could be to chose 
7 o = f B ( ho e ) wo, i.e. rjj = ^ wo- Without this condition, u e 0 belongs only to L p (Fq) for all p £ (2, oo]. 

Therefore, a circulation independent of e and a vorticity with non zero mean value would require to work 
in the Marcinkiewicz space L 2 ’°°(Fq) (weak L 2 space). Even if this space is less classical that L 2 , there is a 
large literature for the well-posedness of the Navier-Stokes equations in fixed domain, because it corresponds 
to relevant initial data for the Euler equations, and also because the self-similar solutions (as the Lamb-Oseen 
vortex) in the whole plane belongs to L 2, °°. In this case, the Cauchy theory is well-known, and Iftimie, Lopes 
Filho and Nussenzveig Lopes managed in [21] to consider the small obstacle problem with non-zero initial 
circulation and initial vorticity with non-zero mean value. For the fluid-solid problem, such a Cauchy theory is 
not yet established. As the optimal decay estimates for the Stokes semigroup are now known in L p — L q [T2] . 
we guess that it would be possible to extend it by interpolation to Marcinkiewicz spaces, and then to prove a 
well-posedness result for the full fluid-solid system. Such an extension would require more work and could be 
interesting, but the main goal of this article is to stay in the standard framework for the fluid solid problem and 
to treat the same question as mm- 

Example 2.5. Another example of initial conditions satisfying (11.1211 can be obtained by truncating a stream 
function associated to a vector field uq defined on R 2 . 

Let us consider uo = V^o 6 L 2 (R 2 ) such that 

divM 0 = 0 in R 2 , curluo £ L 1 D L ? (R 2 ) with q > 1, lim uo(x) = 0. 

| x | —^OO 

We denote by y a smooth cutoff function such that x( x ) = 0 in S(0, 3/2) and x( x ) = 1 in B( 0, 2) c . We consider 
(4,?~o) € R 3 given, then we define 

u E 0 := (i>o(x)x + ( x - X 'j + ro \x-ho\ 

which is divergence free, tending to 0 at infinity, equal to u o far away the solid, and equal to 4 + r o( x ~ h o) 3 * 
in the vicinity of Sq. By local elliptic regularity, we can show that if o € L°°(B(h 0 ,2)) and as tVx (— 7 ^) 
converges weakly to 0 in L 2 , one can check that Uq — L uq in L 2 (R 2 ). 

Actually, as we consider only one solid, we can chose if 0 such that ifo(ho) = 0 and in that case, we can prove 
that the convergence holds strongly in L 2 (R 2 ). 

A last cutoff example comes from the porous medium analysis or the thin obstacle problem (see [TJ US]). 
Instead of truncating the steam function, we truncate directly the vector field uq and we add a correction to 
restore the divergence free condition: 

u o '■= X M x ) + (i - X ^ £ h ° ^j ) (4 + r 0 {x - h 0 ) ± ) + g e (x) 

where g e £ H^(B(ho, 2 e) fl Fq) satisfies 

div g e = - jVx ^ ■ u 0 (x) + ^ • (4 + r 0 {x - ft-o)- 1 ). 

Such a function can be constructed thanks to the Bogovskii operator, and we can prove that 

| 4 e || L 2 < C'e^||u 0 ||t 2 +e|4| +e 2 ko|)- 

See later Proposition 14. 1 1 for details. This implies that Uq uq in L 2 (R 2 ). 


3 Uniform estimates 

Let (uq, 4 i 4) be a family in L 2 (Tq) x I 2 x R verifying ( 11 . 121 ) and such that, up to the extension ( 11 . 101 ) . Uq is 
bounded in L 2 (R 2 ) (see ( 11 . 171 ) 1 . 
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3.1 Change of variables and energy estimate 

As in T2J [5U , we make the change of variables 

v e (t , x) = u £ (t, x — h e (t)), q £ (t , x) = p £ (t, x — h £ (t)), 

and we define 

m = m\t), r e (t) = (o £ y(t). 

The vector fields v £ is the weak solution of a system similar to CE2D (ESI: 


^-+ ([v £ — £ £ ]-V)v £ — diva(v £ ,q £ ) = 0 t > 0, leJJ, 

divu E = 0 t > 0 , x £ Fq, 
lim v e (x) — 0 t > 0 , 


| x | —>-oo 


v E (t, x) = £ £ {t) + r £ (t)x J ~ 

t > 0 , x £ 8 Sq 

m e {£ e )\t) = — f cr(v s ,q £ )n d'y 
■>as~ 

t > 0 , 

J e {r e )'{t) = — f x ± ■ cr(v £ , q £ )n d'y 

JdSS(t) 

t > 0 , 

« E ( 0 ,-) = v £ 0 

in Fq, 

£ £ (0)=t £ o , r £ (0)=r £ o . 



(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 


We set the global density in R 2 : 


p e (x) 


1 x € J-’q, 

p X £ Sq. 


We can define a weak solution 

Definition 3.1. We say that (v £ , £ £ . r e ) is a global weak solution of (EHi-fEU) if for any T > 0, we have 
v £ eL oo (0,T;L 2 (]R 2 ))ni 2 (0,T;77 1 (R 2 )), v £ (t, x) = £ £ (t) + r £ {t)x^ in S £ 0 , 

if it satisfies the weak formulation 

[ p £ v £ ■ (- 7 T— + {[v £ — £ £ ] ■ V)v? eN ) dxds + 2v f f D(v e ) : D(ip e ) dxds = f p £ Vq(x) ■ <p e (0, x) dx, 

J K 2 \ at J J 0 J R 2 Jr2 

for any ip e £ Cf ([0, T); Jt 1 (R 2 )) such that ip e (t,-) £ Vr(Tq). 

One can check that u £ is a weak solution in the sense of Definition o if and only if v e is a weak solution in 
the sense of the above definition. 

We also define the following functional spaces for p £ [l,oo], 

jtp = {v£ £ P (R 2 ) ; divv = 0 in R 2 , D(v) = 0 in S E } , 



with the norm (for p oo) 



For p = oo, the norm £°° 
such that 


is the classical L°° norm. We recall that for any v £ £§, there exists 
v(y)=K+r v y ± (y £ Sq). 


(£ v ,r v ) £ R 3 
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v dy. 


Moreover, one can deduce £ v from v by 

1^0 I fs* 

One can write the system m (Ebd in the following abstract form: 

d t v e + A E v E = divF £ (i> £ ), v e (0) = Vq, 


(3.9) 


where 


V{A E ) := {a € H 2 (M. 2 ) ; diva — 0 in R 2 , D(v ) = 0 in Sq} 
—vAv in J-g, 


A e v := { 2v 


^7 [ D{v)nds+^(f y ± -D( a) 
m J \ J dSn 


ndy\x' L in <Sq . 


(V e 2?(-4 e )), 


F E (v E ) = 


A E := F E A E , 

v s ® (£ v e —v s ) on Tq 
0 on Sq, 


(3.10) 


and where P £ denotes the projector from L P (R 2 ) to and £ v s is defined through (|3.9I) . Note that in the 
definition of A E , D(v)n corresponds to the trace of the restriction of D{v) to the fluid domain. 

The operator —A E is the infinitesimal generator of a semigroup of (S e (t)) t ^o in CA e for p G (l,oo) (see USD- 
Then, Duhamel’s formula gives the following integral formulation of the above equations: 


z (t) = S E (t)v Q+ f S £ (t — s)F E div F E (v E (s)) ds. 

Jo 


(3.11) 


By Sobolev embedding, it is classical to deduce from the weak formulation (see Definition O) that d t v e 
belongs to L 2 (0,T; V' r (Fq)) and that the relation 


f f p E ( d t v e ■ (p s — v E ■ ([a e — £ E ) • V</? £ )) dxds + 2v f f D(v E ) : D(ip s ) dxds = 0 
Jo J r 2 Jo Jr 2 

is satisfied for any ip s G L 2 (0,T; Vr(Fq)). In particular, we can take yf = n e l[o,t], and we remark that 

/ f p E v e ■ ([ v E — £ E \ ■ Vu £ ) dxds = i f f p E (v E — £ E ) ■ V|w £ | 2 dxds = 0 
Jo Jr 2 2 J o Jr 2 

because divu £ = 0 and that [v s — £ e } ■ n|asg = 0. Next, we observe that d t v s G L 2 (0, T; V' r {Fq )) and v s G 
L 2 (0,T;Vr(Fq)) implies that v s is equal for a.e. time to a function C([0,T];£ 2 ) and that 

(!- f p e {x)\v E (t, x )\ 2 dx — 7- f p e (x)\vq(x )\ 2 dx = f f p E (x)v E (t, x) ■ d t v E (t, x) dx for a.e. t G (0, T). 

2 J R2 2 J R2 J Q J R2 

This means that v E satisfies the energy equality 


f p E (x)\v E (t, x)\ 2 dx + 2v f f \D(v s )\ 2 dxds = ^ j p E {x) \vq(x )\ 2 dx for a.e. t G (0, T), 
Jr 2 Jo Jr 2 2 J R2 


hence, we have for any e 


f p £ (f, x) \u E (t, x)\ 2 dx + 2v f f \D{u E )\ 2 dxds = \- j p s 0 (x) \u E 0 (x)\ 2 dx for a.e. t G (0,T). (3-12) 

Jr 2 Jo Jr 2 2 J R 2 


This energy inequality and the hypotheses (11.1411 and (11.1711 imply that 


(■ u E ) e is bounded in L°°(0, T; L 2 (R 2 )) n L 2 (0, T; H\R 2 )). (3.13) 

Let us remark that the energy estimate only implies 

e(h e )'(t), £ 2 {6 e )'{t) bounded in L°°(0, T). 

These estimates do not allow to localize the rigid body and then to use the method developed in Section 01 The 
goal of the sequel of this section is to obtain an estimate of ( h E )' independently of e. 


3.2 Semigroup estimates 

The key for the uniform estimate of i E is the following theorem concerning the Stokes-rigid body semigroup. 

Theorem 3.2. For each q £ (1, oo), the semigroup S E (t) on C q satisfies the following decay estimates: 

• For p £ \q, oo], there exists K\ = K\{p , q) > 0 such that for every Vq £ C q £ : 


II^KHr? < Kit* lq ll u oll,c? for all t > 0. (3.14) 

• For 2 ^ q ^ p < oo, there exists K 2 = K 2 (p, (?) > 0 such that for every F e £ L 9 (R 2 ; M 2x 2 (^)) satisfying 
F e = 0 in 

||S ,£ (t)P e div F e \\ C P < K 2 t~^ + p~^\\F E \\ LQ{R 2 ) for all t> 0. (3.15) 

• For 2 ^ q < 00 , there exists Kg = Ke(q) > 0 such that for every F e £ L 9 (R 2 ; M 2x2 (M.)) satisfying F e = 0 
on Sq: 

Ks E (t)P e div | ^ K(t + ||T ie ||i 9 (R 2 ) for all t> 0. (3.16) 

For e fixed, estimates like (13.1411 - (13.1511 were only established for the Stokes system with the Dirichlet 
boundary condition El Eg. For the fluid solid problem with one rigid disk in R 2 , this result was recently 
obtained by Ervedoza, Hillairet and Lacave in |T3]. The only point to check here is that the constants Ki,K 2 , K? 
are independent of e, which will be easily obtained by a scaling argument. Indeed, as the above estimates are 
optimal i.e. correspond to the decay of the heat solution, they are invariant to the parabolic scaling of the 
Navier-Stokes equations. 

Proof. For e = 1, the statements of the theorem were proved in [12] . see therein Theorem 1.1, Corollaries 3.10 
and 3.11. 

We note that v e (t) := S e (t)vQ satisfies 


dv E 

~dt 


- di vcr(v E ,q E ) = 0 

t > 0, 

x £ 

div v E = 0 

t > 0, 

x £ Ff, 

lim v e (x) = 0 

t > 0, 


|cc|—>-oo 



= f(t) + r E (t)x ± 

t > 0, 

x £ DSq 


m e {l E )'{t) = — f cr(v e ,q E )n dj t > 0, 
J e (r E )'(t) = — f x 1 ' ■ cr(v E , q e )n d'y t > 0, 

J dSZ (t) 


r (0,-) = ' 


7T< 
J ( 


0 5 


£ e (0 )=£§, r E (0)=r E o . 
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Setting 


(3.17) 


v(t, x) := v e (s 2 t,sx), q(t, x) := eq E {e 2 t, ex), £(t) := £ e (e 2 t), r(t) := er E (e 2 t), 
standard calculation gives that 


dv 

dt 


-' 7 — — div a(v, q) = 0 t > 0 , Jq 1 , 


div v = 0 t > 0 , x £ Fq 


0) 


lim v(x) = 0 t > 0, 

| X | —^OO 

v(t, x) = £{t) + r(t)x ± t > 0, x G 8 Sq, 

mn}£'(t) = — / a(v,q)n dj t> 0, 

JdSl 

J l r'(t) = — x J ~ ■ a(v,q)n dj t> 0, 

JdSlft) 

u(0,-)=W in-Fo 1 , 

*( 0 ) = 4 , r( 0 ) = r 0 , 


where 


v 0 (x) := Vq(ex), £ 0 := £ e 0 , r 0 := er e 0 . 
This means that v{t) = u 1 (t) = S' 1 (t)r >0 and thus that 

||u(t)|| £ p < ||uo|| £ 9 for all t > 0. 

Using (|3.17[1 - (|3.18[) . this estimate is equivalent to 

\\v e (t)\\c* < K it*~«\\vo\\d for all t> 0. 

Relations (13.151) and (13.161) can be done similarly. In that case, we also set 

F(x) := -F e (ex) 

£ 

and we show that if v E (t) = S e (t) P e divU 6 , then v defined by (13.171) satisfies 

v (t) = S\t) P 1 div F. 


(3.18) 


□ 


3.3 Uniform estimate on the solid velocity 

We first show that there exists Ao > 0 such that if ||uq||,c 2 ^ Ao, then there exists a unique 

G C°([0,T];£ 2 ) nC 3 ° / 8 ([0, £];£«) with £ E := £ v e G C° / 2 ([0, T]; R 2 ) (3.19) 

satisfying (13.111) (that is a mild solution). Here we have denoted for any Banach space X by C°([0, T]; Al) the 
Banach space of functions / such that 1 1 —>■ t a f(t ) are continuous from [0,T] in X. The norm associated is 

ll/llc°([ 0 ,T];X) : = SU P ^Wfi^Wx- 

te[o,T] 
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Proposition 3.3. There exist Ao,/ro > 0 independent of £ such that the following holds for any T > 0. If 
Vq € C 2 satisfies 

IKIUj < Ao (3.20) 

then there exists a unique v e satisfying (13.111) and such that 


IKIIc°([o,t];£ 2)> ||« £ || c o /8 ([ 0 ,T];£f). IKIIc° /2 ([o,t] ; rU are bounded by 


Moreover there exists a constant C > 0 independent of £ such that, if vo a > v ob 
satisfying (13.201) . then 


IK - « 6 llc°([ 0 ,T];£*) < C\\v e 0a - ug b || £ 2 . 


G C 2 are two initial conditions 


(3.21) 


Proof. As we have proved in the previous theorem that the constant in the semigroup estimates are independent 
of £, it is then enough to follow the fixed point argument in [T2J pp. 364-371]. For completeness, let us write 
here the details. 


Let us introduce the space 

X s :={v s eC%[0,T)-,jr 2 e )nC° 3/s ([0,T]-Cl) with t v * eC 1 ° /2 ([0,r];M 2 )} 

endowed with the norm 


u ■= 




o v L 


Let us also define the map 
defined by 


X s , 


Z E (v E )(t) = S E {t)v E 0 + [ S E {t-s)F E dWF E (v E )(s)ds, 

Jo 

where F e is defined by (13.101) . One can define 

${v E ,w E )(t) = f S E {t- s)F E dWG E (v £ ,w E )(s)ds, 

Jo 

where 

r e ( ,, s ® _ wE ) 011 T o 

1 ’ ’ \ 0 on 5g. 

We deduce from (13.151) that 

ti\\d>{v E ,w E )(t)\\ci f (t-s)" 5/8 (||w £ (s)®w e (s)|| £ 4 + |^ e (s)|||n e (s)|| £ 4) ds. 

Jo 

Using Holder’s inequalities, we obtain from the above inequality that 

ti||$(u e ,u; e )(t)||£| ^tiK 2 ( 8,4) f (t - s)~ 5/s (s~i |Kllc° /8 £f s~i ||w e || c ° /8 .ci 

J 0 

+ s-5|^ E | c o /2 ||^||^i £2 (s-|) 2/3 ||u e ||^ 3 g£i ) ds 
<2if 2 (8,4)B(| J ^||t, E |U.||«, e |U., (3.22) 

where H(-,-) is the Beta function: 

B(a,/3):= f (1 — r)~ a T ~ 15 dr. 

Jo 
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Similarly, 

Finally, 


\\$(v s ,w s )(t)\\ cl ^2K 2 (2,2)B 

L 


i i 

2 ’ 2 


,W £ )(t) 


= / 4 


F 


S £ (t—s) P e div G e (v e ,tt) e )(s) 


and from ()3.16j) . we deduce 

/ t 

Ke{4)(t- s) _ (5 + 3)||G e (u e ,M; e )(s)|| L 4 (R 2 ) ds. 


With the same estimates as in (13.2211 . we obtain 


3 3 


t 2 \h( v e,w‘)(t)\ < 2K t (A)B IKIMKII**- 

Gathering (13.2211 . (13.2311 and (13.2411 yields 

ll*(«‘> e )IU- < ColKIMKIU-, 

where 


5 3 


8’ 4 


Co = 2 K 2 (8,4)B -, - + K 2 (2 1 2)B -, - + A',(4)B 


1 1 


2 ’ 2 


3 3 
4’4 


We assume (13.2011 for Ao that we fix below and we apply (13.1411 in order to obtain 

||S e i;§||*. <CrA 0 , 

where 

Ci=iG(8,2) + iG(2,2) + Arr(oo, 2). 

Note that we have used in (13.2611 the relation 

which is a consequence of (ESI. 

Relations (13.2511 and (13.2611 imply that the mapping Z e is well-defined and that 

\\Z s {v*)\\ X e <CiA 0 + q,||w e |&.. 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


Let us set 


R := ——— and Ao = min 
4Gn 


f— 

V 2 C 1 


,r)- 


(3.27) 


Then the closed ball Bx^(0,R) of X e is invariant by Z e and if v e ,w E G Bx*(0, R), 


II Z e (v e ) - Z e (w e )\\ X e = ||$(w e ,w £ - w E ) + $(u e - w e ,w e )\\ X e < 2C 0 R\\v e - w e \\ x * = ^IK “ w e |U=- 


Using the Banach fixed point, we deduce the existence and uniqueness results. Moreover, we have p,o = R = 
l/(4Go), which is independent of e. This strategy comes from |22] and [23]. It is originally done through an 
iterative method, but it was adapted as a fixed point argument in [2S] (see also p?2]'). 

Sensitivity of v e to the initial data. Assume VQ a ,Vg b G satisfy (13.201) with Ao defined in (13.271) . Then, 

\\ v l~ v b\\x* < II^Ko ~ v 0b)\\xe + \\^K,v £ a - vl)\\ xe + \Mv e a - vl,vf,)\\ xe 

< Cill^a ^ + Co||< - vf\\xe (|K|U= + \\vt\\xe) 

< Cl||Uoa - v 0 b\\cl + 2 C 0 M 0 IK - vl\\ X e. 

We conclude that 

IK - v t\U‘ < 2 Ci||^ a - vfoWc,. 

□ 
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One can show that a mild solution in the above sense is also a weak solution (see [12]). For sake of 
completeness, we give the proof of this result here. 

Lemma 3.4. Assume V s satisfies (13.19D - (13.20D and (13.111) . Then v e is the weak solution of (13.II) - (13.81) in the 
sense of Definition \3.1\ 

Proof. Let us consider a sequence (vf in ) n with values in V^A 6 ) 1 / 2 ) such that 

v e on -> Vo in C 2 e 

and such that Vq u satisfies (|3.20l) . For all n, it is proved in [21] that there exists a unique strong solution 

< G H 1 (0, T; Cl) O C([0, T]]V{{A e y/ 2 )) n L 2 (0, T; V{A e )) 
of (13.1D - (13.8D and it satisfies (13.111) since in this case 

divF e «)£L 2 (0,T;L 2 (R 2 )). 

It is also proved in [21] that (u®) converges towards the weak solution of (I3.1D - (I3.8D in the sense of Definition l3.ll 
From (13.211) . we also have that (u£) converges towards the mild solution of Proposition 13. 31 in L°°(0, T; £ 2 ). 
Consequently, the mild solution v e , that satisfies (13.19D - (]3.20D . is the weak solution associated to Vq. 

□ 


4 Proof of Theorem 11.2 

This section is dedicated to the proof of Theorem 11.21 As recalled in the introduction, [21] established that for 
any e > 0, there exists a unique weak solution (u e , h s ,0 e ) to problem m cu) in the sense of Definition 11.11 
Let us fix T > 0. 

4.1 First convergences 

Thanks to the energy estimate (13.131) . we can extract a subsequence such that 

u E ^u in L oo (0,T;L 2 (R 2 ))nL 2 (0,T;ll 1 (R 2 )). (4.1) 

By abuse of notation, we continue to write u e the subsequence. 

If the initial data satisfies the smallness condition (11.161) with Ao given in Proposition 13.31 we deduce from 
Section 13.31 that 

|(/0'(t)|<^ (t > o), 

where fj ,o is independent of e. As a consequence, 

\h e (t)\ ^2[i 0 Vr (t G [0,T]). 

We fix q G (1,2), thus ( h e ) is bounded in W 1 ' q [f) 1 T;M 2 ) and we have, up to a subsequence, 

h e —> h uniformly in [0, T], (4-2) 

with 

h€W^ q {Q,T). (4.3) 
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4.2 Modified test functions 


The key to treat shrinking obstacles problem is to approximate test functions 
in the perforated domain. 

in R 2 by admissible test functions 

Proposition 4.1. Let T > 0, p £ C%°([0,T) x R 2 ) with divy> = 0. We consider q £ (1,2) as in Section \4~1\ 
For any rj > 0, there exists p v £ W c 1,9 ([0,T); H 1 (M. 2 )) satisfying 

divp r > = 0 in [0, T) x R 2 , 

(4.4) 

<P V = 0 t £ [0, T), xeB(h(t),?), 

(4.5) 

p^^p L°°(0, T; JT 1 (M 2 )), 

(4.6) 

d t p v ->■ d t p L q (0, T ; L 2 (R 2 )). 

(4.7) 


Proof. We introduce a cutoff function y G C°°(R 2 , [0,1]) such that \ = 1 in B{ 0, l) c and y = 0 in B( 0,1/2). 
Let us denote the annulus B( 0,1) \ B( 0,1/2) by A. 

We remark that the function dp 17 : (t, y) H>• </?(t, yy + h{t))'S/ x(y) belongs to W 1,<? (0, T; L 2 (A)) and verifies for 
any t 


P r ’{t,y)dy= di v(ip{t,riy + h(t))x(y))dy= p(t,yy + h{t)) ■ n{y) ds 

1 J A v 7 JdB(0,l) 

div r/y + h{t))j dy = 0, 


'5(0,1) 


where we have used twice that <p is divergence free, that y = 1 on dB{ 0,1) and that y = 0 on dB{ 0,1/2). With 
these properties, it is known by m Theorem III.3.1] (and Exercice III.3.6) that there exists C depending only 
on A such that the problem 

div g v = (p n , g v £ W 1,q (0,T; Hq(A)) 

has a solution such that 


llff ?? l|L-(0,r;Bi(A)) < C'||^ ?, ||l°°(0,T;L 2 (A))> 

ll^t5 T, ||L < '(0,T;B 1 (A)) ^ C||5i^ ?, ||L9(0,T;L 2 (A))- 


Extending g v by zero in the exterior of A, we define 


p v (t, x ) = p(t, x)x( x - g v (t, x) 


where 


9 




We easily verify the divergence free condition m- Moreover, with a change of variables, we also note that 
1 


9 


5 


L“(0,T;L 2 (R 2 )) 


V.g r 


L°°(0,T;L 2 (R 2 )) 


^ C'||^ ?7 ||L«’(0,T;i 2 (A)) ^ C1MU“((0 ,T)xR 2 ) 


(4.8) 


so we check that 


~\\ l P V - 7 , IU°°(0,T;L 2 (R 2 )) + ||- Vv5||ioo( 0 ,T;L 2 (R 2 )) < C'|M|w 1 .°“((0,T)xR 2 ) 

which gives directly that p v converges to p strongly in L°°(0,T; L 2 (R 2 )) and weak-* in L°°(0, T; H 1 (M. 2 )). By 
uniqueness of the limit, we do not need to extract a subsequence and we get the weak limit (|4.6I) . 
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Now we compute 


dt^i^x) - d t p{t,x) =d t ip(t,x)(x( X - l) 


vfox) wu\ tx~7 \( x ~K t ) 


V 

d t g v (t, X + h'(t) ■ Vg ri (t,x). 


'{t) ‘ (Vx)(- 


It is obvious that the first right hand side term converges to zero strongly in L°°(0, T; L 2 (R 2 )). It is also an 
easy computation to check that the second term is bounded in L 9 (0, T; L 2 (R 2 )) and tends to zero strongly in 
L 9 (0, T; L P (R 2 )) for p £ [1,2). Hence, it converges weakly to zero in L 9 (0, T; L 2 (R 2 )). From (Id.31) and (Id.81) . we 
know that the last term is bounded in L q (0 , T; L 2 (R 2 )), and as it converges to zero in T ) x R 2 ), we infer 

that it converges also weakly to zero in L 9 (0, T; T 2 (R 2 )). Finally, we note that 


1 

V 



x — h(t) 
V 


L°°( 0 ,T-L 2 (R 2 )) 


^ C||<9 t <p||Loo(( 0jT ) xR 2), 


hence the third right hand side term tends to zero strongly in L°°(0, T; L 2 (R 2 )). This gives (Id.71) . 

Due to the support of x and g v , it is clear that ip v = 0 on B(h(t), §). This ends the proof. □ 

Remark 4.2. An important consequence is the approximation of any test function. Let T > 0 and p £ 
C^°([0,T) x R 2 ) with divtp = 0. Then, we have constructed a family [p r, ) r i of divergence free test func¬ 
tions which tends to p in the sense of (ld.6D - (Id.71) . Moreover, for any rj > 0 fixed, we put together the strong 
convergence of h £ (O with the support of p v (USD to deduce the existence of e v > 0 such that 


p v = 0 fortEp,! 1 ), x £ S e (t), £^£n- 

This implies that p v is an admissible test function for the fluid-solid problem (see Definition II.ID . 

Remark 4.3. In the proof of the above proposition, we note that H 1 is the critical space in dimension two: 
Xijj) ~ 1 tends to zero strongly in W 1,p for any p £ [1, 2), is bounded in H 1 (then tends weakly to zero), and 
goes to infinity in W 1,p for p > 2. This explains why the standard framework for shrinking obstacles problems 
is H 1 (see, e.g., (HE El [201 [2D El)- Nevertheless, as we need an estimate of the solid velocities, it is natural to 
look for a C° estimate of the velocity, hence a H s estimate for s > 1. Unfortunately, the H 2 analysis developed 
in mm fails (see 0 ). 

In dimension three, the critical space for the cutoff argument is IF 1,3 which is again not embedded in C°. 


4.3 Passing to the limit in the Navier-Stokes equations 

The first step is to pass to the limit e —> 0 for p fixed. 

Theorem 4.4. Let T > 0 and let p £ C^°([0, T ) x R 2 ) with div ip = 0. We consider the family (<p I? )r;>o obtained 
in Proposition o Then, for any rj > 0, the limit u of u e (see 63; verifies 



^ + (». v V 


dx + 2v f f D{u) : D(ip v ) dx 
Jo Jr 2 



p v (0, x) dx. 


Proof. Let g > 0 be fixed. From (Id.3D and Sobolev embeddings, we know that t h(t) is continuous on [0, T] 
and then uniformly continuous. Hence, there exists a uniform subdivision to = 0 < t± < ... < 4 m+i = T such 
that for any t £ (tj, ij+i), 

I h(t) - h(tj )| < |. 

From (ld.5D in Proposition ld.il we deduce 

= 0 in {tj,t j+ 1 ) x |). 


15 











Putting together this relation with (14. 21) . there exist open relatively compact sets Oj and > 0 such that for 
all e < i v 

M 

S e (t) (~1 Oj = 0 for all t £ (tj,tj + 1 ) and supp^ C ^^(tj,tj+ 1 ) x Oj. (4.9) 

j=o 


For any j = 0,..., M, we write the Helmholtz-Weyl decomposition 


u E = p OjU 6 + Vg e , 

where Pe^ is the Leray projection on H(Oj) (see the introduction for the definition of 71(0)). This projection 
is orthogonal in L 2 and by a standard estimate on the Laplace problem with Neumann boundary condition, 
there exists a constant Co 3 - > 0 such that 

l|Po i M e ||i 2 ( C)i ) ^ \\u e \\ L 2( 0 .) and \\PojU e \\m(Oj) < C 0j \\u e \\ H i( 0 .y 

Thus, by (13.131) . 

(P 0j u e ,S7q e ) e is bounded in L°°(0, T; L 2 (Oj)) D L 2 (0, T; H x (Oj)). 

In particular, 

P 0j u e ^P 0j u in L°°(0, T ; L 2 (Oj)) D L 2 (0, T ; H 1 {O j )) ) (4.10) 

Vq e ^ Wq = u-p 0j u in L°°(0,T-,L 2 (O j ))nL 2 (0,T;H 1 (O j )). (4.11) 


Now we derive a time estimate for P OjU £ in order to get a strong convergence. For any divergence free test 
function ip £ C%°((tj,tj+ 1 ) x Oj), we have by (14.91) that ip(t,-) £ Vij(.F e (t)) (see (11.111) '). hence (11.131) gives 

r t j +1 r 

(a t P 03 .u £ ,^) L 2 ((t . it . +l);V(0 . ) , )ji 2 ((t . it . +l);V(0 . )) =- / / P 0j u e ■ d t ipdxdt 

Jti Jo< 


tj J Oj 
rb+i r 

Jtj JOj 
r t 5 +1 


• dtip dxdt 


= f f u e ■ (u e ■ S7)ip dxdt — 2v f ( D(u e ) : D(ip) dxdt. 

J t j J Q j J tj J G j 

Thus, bv using (13.131) and the internolation ineaualitv II f llr^iim ^ II ||V/||^ 2 2 (r2) , we get 


(0 t Pc. 3 .U £ ,^)L2((t j ,t j+1 );V(O j )'),i 2 ((L>b+i);V(Od) 

< ll u 1lL4(( tj> t J . +1 );L4(od)ll^lli 2 ((b.b+i);V(C'd) + \\ Du£ \\L^(tj,tj +1 y,L^(o j ))\\i’\\LH(tj,t j+1 )-v(Oj)) 

< C'IIVi|L 2 ((t 3 -,t 3+1 );V(O j ))- 

Consequently, (9 t P OjU e ) is bounded in L 2 ((tj, tj+i); V(Oj)'), and the Aubin-Lions lemma in H 1 D 71(0j) °->- 
L 4 nn{Oj) V(Oj ) allows us to extract a subsequence such that 

POjU e ->• Po,. u strongly in L 2 ((tj,t j+1 )-, L 4 (Oj)). (4-12) 

Actually, by the uniqueness of the limit, we do not need to extract a subsequence in (14.121) . 

These convergences are enough to pass to the limit in the Navier-Stokes equations. Indeed, for any e £ (0, s v ], 
we know from (TOl) that ip v is an admissible test function, and (11.131) reads 


Jo Jv at £- 


j =o 


(u e ® u e ) : Vip 11 dxdt + 2u [ f D(u e ) : D(p v ) dxdt 

Jo Jr 2 

= f Uq(x) ■ p v (0, x) dx. 

Jr 2 
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Using the weak limits (14.111 and (11.171) . we easily pass to the limit in the linear term 


u" • dxdt + 2v [ [ D(u e ) : D(<p v ) dxdt 


1 0 J R 2 


dt 


o Jr 2 


ff 

Jo Jr 2 


r T r d(pV rT 

u ■ dxdt+2v 


dt 


n 

Jo Jr 2 


D{u ) : D(ip") dxdt 


and 


/ Uq(x) ■ <p v (0,x) dx — > j uq{x) ■ <p v (0,x) dx. 

Jr 2 Jr 2 

For the non-linear term, we decompose in (tj,tj+ 1 ) x Oj as 

u £ ® u £ = Poy« £ ® u e + X7q £ ® P OjU £ + S7q £ ® Vg £ . 

Let us note that for any harmonic q (i.e. A q = 0), we have the following relation: 

/ (Vg ® Vg) : Vip" = - f div(Vg ® Vg) • ip v = - f ( iv|Vg| 2 • cp" + AgVg ■ ip" J = 0, (4.13) 

J Oj J Oj J Oj \ / 

because ip" is divergence free and compactly supported in Oj. From (14.10f> - (|4.12|> . we have 
/ (P OjU e ® w £ ) '■ Vtp" dxdt —> / (P OjU® u) :Vp" dxdt , 

■ (tj 4j+i)xOj J(tj,tj+ 1)x Oj 


and 


I (Vg £ ® ¥ojU £ ) : 'S/p" dxdt — > I (Vg ® Po 3 .u) : S/p" dxdt. 

i)xOj J (tj ,tj+±)xOj 


Gathering the two above convergences and (14.131) applied to q e and to g, we conclude 

/ (u £ <8) u £ ) : V p" dxdt —»■ / (u ® it) : Vg? 17 dxdt. 

(tj ,tjj-i ) x Oj J {tj ,tj+±)xOj 

This ends the proof of Theorem 14.41 

To end the proof of Theorem ll.2l it is sufficient to pass to the limit ij —> 0, thanks to Proposition l4.il 


□ 


Proof of Theorem Al.Sl Let T > 0 and p £ (^([OjT) x R 2 ) with divy> = 0 fixed, then we consider {p") v ^rn 
which approximate p (see Proposition 14.11) and u a weak limit of u e (see (14.11) 1. Theorem 14.41 states that the 
limit u verifies for any iq 


■ff 

Jo Jr 2 


(dip" 

— 


■ SZ)p"'\ dx + 2v f f D(u) : D{p") dx = f uo{x) ■ p"(0, x) dx. 
J Jo Jr 2 Jr 2 


/o Jr 2 V ^ 

As u belongs to L°°(0, T; L 2 (R 2 )) D L 2 (0,T; H 1 (R 2 )), we deduce from the convergences (14.6D - (I4.7I) of p" to ip 
that 

— f f u ■ ( -ff- + (u ■ VV) dx + 2v f f D{u) : D(ip) dx = f uq(x) ■ ip( 0, x) dx. 

Jo Jr 2 \vt J Jo Jr 2 Jr 2 

By density, this equality is also true for any ip £ ([0, T); V(R 2 )). Noting that 

J D{u) : D[ip) = — J V« : S/p + - J divrtdivg? = — J V« : S7ip, 


we conclude that u is a weak solution to the Navier-Stokes equations in R 2 associated to uq. By uniqueness of 
such a solution, we note that the weak convergence (ED holds for all sequence (e n ) converging to 0: 

u £n ^u in L oo (0,T;L 2 (R 2 ))nL 2 (0,T;iL 1 (R 2 )) 

as n —> oo, without extracting a subsequence. 

This ends the proof of Theorem 1 1.21 □ 
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5 Final remarks and comments 

5.1 The three dimensional case 

The main obstruction to the generalization to the three dimensional case is that the optimal decay estimates of 
the Stokes semigroup are not established for all time t > 0. 

For a solid with any shape moving in a three dimensional viscous fluid, Geissert, Gotze and Hieber m 
Theorem 4.1] show the maximal regularity for the Stokes semigroup, locally in time. By an extension operator 
and Sobolev’s embedding, we can deduce from their result the optimal L p — L q estimates for some p , q (see 
Em Proposition 3.1]). Even if we have to check that the p, q reached are enough to perform the fixed point 
argument (as in Proposition l3.3D . the problem is that these estimates are local in time, i.e. valid up to a time T, 
with the constants depending on T. Therefore, after the scaling argument (see the proof of Theorem l3.2D . these 
estimates are independent of £ only up to a time T e = eT. Indeed, by the scaling property of the Navier-Stokes 
equations, the small obstacle problem is similar to the long-time behavior. 

In dimension two, the global optimal L p — L q decay estimates for the Stokes semigroup were obtained 
by Ervedoza, Hillairet and Lacave in |T2] when the solid is a disk in the whole plane. We guess that their 
analysis can be adapted in the exterior of a ball in dimension three, but it would require a considerable work, 
decomposing the Stokes equations on spherical coordinate system instead to polar decomposition, exhibiting 
the “good unknown” (see [HI Proposition 2.3] for the two dimensional case), and adapting the elliptic lemmas. 

Finally, one should also adapt the fixed point argument performed in Proposition 13.31 which should not be 
too difficult, because the original proof of Kato [22] holds for any dimension n ^ 2. 


5.2 Extension to massive pointwise particles 

If some solids 

Sf ]q := htfi + £iSi,o, 

shrink to massive pointwise particle, i.e. 

ml = m\ > 0 and Jf = e 2 j} > 0, 


(5.1) 

(5.2) 


then the energy estimate gives directly the uniform estimates of the solid velocities. Therefore, we do not need 
here the analysis developed in Sections 13.1113.31 and we can prove the following result for any shape, and with 
several solids. 


Theorem 5.1. Let N rigid bodies Sf(t) of shape Si, q (15.11) (with Si,o smooth simply-connected compact subset 
ofM. 2 , with nonempty interior and where the center of mass of Si,o is 0) inside a bounded domain LI, and where 
the positions hi ,o £ LI and the size e are chosen such that 


s to n S j,o = 0 (* ^ J) and Sf 0 C LI. 


Let Uq £ Vr(Tq) such that 


Uq uq m 


L 2 (Ll ), 


Kf.o I 


c , e|r? 0 | < C, Vie {!,..., TV}. 


(5.3) 


There exists a global weak solution (u e , hi , 9() see I ITtM (see also my, ®). Then there exists T > 0 such that 


u e ^u in i°°(0, T; L 2 (Ll)) DL 2 (0,T; Hq(LI)) 
where u is the weak solution of the Navier-Stokes equations associated to uq in LI. 

Remark 5.2. To adapt Proposition 14.II we need the existence of a positive distance between the rigid bodies, 
independent of e. Therefore, denoting by hi the limit of hi, the time T in the above theorem corresponds to a 
time such that 

\hi{t) - hj(t)\ > 0 {i^j, ie[0,T]). 

The existence of such a T is ensured by (15.311 . 
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Namely, we can prove the following: 


Proposition 5.3. Let T > 0, ip £ C“([0,T) x f2) with div</? = 0 and consider ip > 0 such that 

| hi(t) — hj(t )| ^ 2r/i for all t £ [0, T\ and i / j 


and 

dist(supp ip(t, •), dtt) ^ 2r)i for all t £ [0, T], 

For any 77 ^ r)\ there exists gf 1 £ W c 1 ,oo ([0, T); Hq (LI)) satisfying 

div gp = 0 in [0, T) x SI, 

^ = 0 t £ (0, T), x £ B (hi(t), |) , 

d t ^^d t g> L°°(0, T ; L 2 (fl)). 

In contrast, we do not need a positive distance between the bodies and the exterior boundary dLl. In 
particular, in the case of a single rigid body (i.e. N = 1), we can take T arbitrary large. 

The rest of the proof can be done following Section 14.31 
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